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By a novel exact diagonalization technique we show that bound pairs propagate 
between repulsive Hubbard clusters in a superconducting fashion. The size of the 
matrices that must be handled depends on the number of fermion configurations per 
spin, which is of the order of the square root of the overall size of the Hilbert space. 
We use CUO4 units connected by weak O-O links to model interplanar coupling and 
c-axis superconductivity in Cuprates. The numerical evidence on CU2O8 and CU3O12 
prompts a new analytic scheme describing the propagation of bound pairs and also the 
superconducting flux quantization in a 3-d geometry. 



Evidence for pairing in the repulsive Hubbard and related models has been reported by several authors. 
Analytic approaches [1] [2], even at strong coupling [3], generalized conserving approximation theories like 
FLEX [4], as well as Quantum Monte Carlo Studies on supercells [5] lead to this conclusion. However, we 
want more evidence about the real nature of the pairing interactions. Su and coworkers [6] have reported 
that in narrow bands one- and two-dimensional Hubbard models no kind of superconducting long-range order 
holds at any non-zero temperature. Here we wish to explore the possibility that the Hubbard model can 
show superconductivity in the ground state when interplanar coupling is allowed. Since one cannot master the 
problem with an infinite stack of infinite planes, some economy is needed. However in high-T c superconductors 
the coherent length is ~ a few lattice constants, and Cu-O planes can be approximately represented by 
clusters that are large enough to host a bound pair. The interplanar hopping does not dissolve pairs and 
superconducting flux quantization is their clear signature. The magnetic properties of attractive Hubbard 
models have been studied by Canright and Girvin [7]; here we propose a gedankenexperiment very much in 
the spirit of Little and Parks [8] , in the repulsive case. 

The repulsive Hubbard Hamiltonian of fully symmetric clusters C has two-body singlet eigenstates without 
double occupation [9] [10] [11] [12] called W = pairs. The presence of such solutions at the highest occupied 
level of the non-interacting (Hubbard U — > 0) system is necessary to allow Ac(N) < where Ac(N) = 
E { C 0) (N) + E ( °\N - 2) - 2E { C 0) (N - 1), and (N) is the interacting ground state energy of the cluster C 
with N fermions. By means of a non-perturbative canonical transformation [2] [13], it can also be shown 
that Ae(iV) < is due to an attractive pairing effective interaction and at weak coupling |Ac(A r )| is just the 
binding energy of the pair. 

CUO4 is the smallest cluster which fully preserves the point-symmetry of the Copper-Oxide planes of high- 
Tc materials. We have already described W = pairing in great detail as a function of the one-body and 
interaction parameters on all sites; the study was extended to larger clusters too [1] [11]. W = bound pairs 
in the CUO4 cluster are found to exist in the physical region of the parameter space. However, since it is the 
symmetry that produces the pairing force we use the simplest working model to study bound pair propagation. 
Here, in order to simplify the analytical formulas, we neglect the 0-0 hopping term and also any distinction 
between Cu and O sites (except geometry, of course). The only nonvanishing hopping matrix elements are 
those between an Oxygen site and the central Copper site; they are all equal to t. For the sake of simplicity, 
we parametrize the Hubbard model in such a way that actually everything depends only on the ratio U/t; 
the important thing is that in this way we still have access to the part of the parameter space where pairing 
occours [9]. Thus, we consider the Hubbard Hamiltonian 

#Cu0 4 =tJ2(dip™+Pld°) + U(Y,^ ) ^ > +nf ) n{ d) ) (1) 

ia i 

where p\ a and pi a are the creation and annihilation operators onto the Oxygen i = 1, .., 4 with spin a = j, {, a% 
and dfj are the creation and annihilation operators onto the Copper site, while nfj = p\ a Pia and ha = d\d„ 
are the corresponding number operators. -Hcu0 4 is invariant under the permutation Group S[4], which has the 
irreducible representations (irreps) Ai (total-symmetric), B2 (total-antisymmetric), £ (self-dual), 71 and its 
dual 7-2, of dimensions 1, 1, 2, 3 and 3, respectively. The ground state of Cu04[2] (i.e. CUO4 with 2 fermions) 
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belongs to x Ai and that of CuC>4[4] is in £; both are singlets, as the notation implies. The ground state of 
Cu0 4 [3] is a 2 Ti doublet. A Cu o 4 (4) < for this model when < U < 34.77 t, as shown in Fig.(l). 



A Cu0 (4) 




FIG. 1.: Trend of Ac u o 4 (4) in t units versus log[U/t] 



Thus, we introduce a graph A with CUO4 units as nodes. The total Hamiltonian is 

Htot ~ Ho + H T . 

with 



(2) 



(3) 



where itT is the creation operator onto the Oxygen i = 1, .., 4 of the a-th cell and so on. Hence, the point 
symmetry Group of Ho is 5[4]' A ', with |A| the number of nodes. There are many different ways to model an 
inter-planar hopping. Nevertheless, to preserve the symmetry that produces the Acu0 4 (4) < property, H T 
must be invariant under the S[4] diagonal subgroup of S[4] |A| . In the following we shall consider a hopping 
term allowing to a particle in the i-th Oxygen site of the a-th unit to move towards the i-th Oxygen site of 
the /3-th unit with hopping integral r a 0: 



T P/3,ia 



+ h.C.] 



(4) 



a,0GA 



For N = 2|A| and T a p = 0, the unique ground state consists of 2 fermions in each CUO4 unit. This paper 
is devoted to the inter-planar hopping produced by small r Q /3 <C |Ac u o 4 (4)| with a total number of particles 
N = 2| A| + 2p; p represents the number of added pairs. When U /t is such that Ac u o 4 (4) < 0, each pair prefers 
to lie on a single Cu04 and for N = 2|A| + 2p the unperturbed ground state is 2 p x(' A ') times degenerate 
(since 1 £ has dimension 2). 

By this sort of models one can study the interaction of several fermion pairs in the same system. The 
simplest topologically non-trivial graph is the ring, with a set A = {1, 2, . . . ,|A|} and 



( t if /3 = a + l, 
T a p = < t* if p = a — 1, 
I otherwise . 



(5) 



where (f> is the magnetic flux concatenated by the ring and (j>o = — . In the absence of magnetic field, r will 
be taken to be real. 

Note that for p=0 the concentration (number of holes per atom) is 2/5 = 0.4; this is somewhat more than 
half-filling (1/3 w 0.33) but still reasonable. We are using CUO4 as the unit just for the sake of simplicity, but 
the W = mechanism produces bound pairs at different fillings for larger clusters [12] and the full plane [1] 
[2] too. By replacing Cu04 by larger units one can model other ranges of the hole concentration. 

We exactly diagonalize the |A| — 2 and |A| = 3 ring Hamiltonian; to this end we introduce the Spin- 
Disentangled technique. We let Mj + Mj = N where M a is the number of particles of spin a; {\4>aa)} is a real 
orthonormal basis, that is, each vector is a homogeneous polynomial in the p^ and of degree M a acting on 
the vacuum. We write the ground state wave function in the form 



a/3 



L a p\4>a1;) <8> \<j>pi) 



(6) 



which shows how the f and J. configurations are entangled. The electrons of one spin are treated as the "bath" 
for those of the opposite spin: this form also enters the proof of a famous theorem by Lieb [14]. In Eq.(6) L a( g 
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is a mj x mj rectangular matrix with W CT =(^ )• We let K a denote the kinetic energy rrv x m a square matrix 

of Htot in the basis {|</w)} 5 and ArJ the spin-cr occupation number matrix at site s in the same basis (iV s 
is a symmetric matrix since the |<^w)'s are real). Then, L is acted upon by the Hamiltonian i/tot according 
to the rule 

Hu*[L] = [KtL + LK{\ + UJ2 N ^LN ( S 1) . (7) 

S 

In particular for Af T = Mj. (S z = sector) it holds K r = K l and /Vj T) = iV s (i) . Thus, the action of H is 
obtained in a spin-disentangled way. In the S z — sector for |A|=3 the size of the problem is 1863225 and 
the storage of the Hamiltonian matrix requires much space; by this device, we can work with matrices whose 
dimensions is the square root of those of the Hilbert space: 1365 x 1365 matrices solve the 1863225 x 1863225 
problem, and are not even required to be sparse. We believe that this approach will be generally useful for 
the many-fermion problem. Since we are mainly interested in getting the low-lying part of the spectrum as 
fast as possible we opted for the Lanczos method, taking advantage from our knowledge of the S[4] irrep of 
the r = ground state; the scalar product is given by (*!|* 2 > = Tr{L\L 2 ). In this way, the Hamiltonian 
matrix takes the tri-diagonal form; however a numerical instability sets in well before convergence is achieved 
if one uses chains longer than a few tens of sites. Therefore we use repeated two-site chains alternated with 
moderate-size ones. 

The two-Cu04 ring (14,400 configurations in the S z — sector) is readily solved by a Mathematica code on 
the personal computer; however, this cluster is not adequate for studying the quantization (superconducting 
or otherwise) of a magnetic flux by the bound pair. The reason is that the two units are each at the left 
and at the right of each other; any vector potential perpendicular to the CuCU's can always be gauged away. 
However, we have verified that the ground state energy with 6 holes is ^c°u 2 o s (6) — ^Cu0 4 (4) + ^Cu0 4 (2) f° r 
t — and it receives a negative correction oc r 2 /|Ac u o 4 (4)| for small r, which is consistent with the presence 
of a bound pair. 




FIG. 2.: Numerical results for CU3O12 with r = 0.001. Lowest-energy eigenvalues labeled by their interplanar quasi- 
momentum are shown versus flux <j>. The pattern is periodic (a flux quantum can be gauged away), (a) U = 0. A 
paramagnetic current is excited by the field and the system is utterly normal. (6) (7 = 5. The ground state shows a 
clear superconducting pattern, with a minimum at <j> = 4p- . All energies are in t units. 

The three-Cu04 ring behaves similarly, but can also concatenate a flux. In Figure (2. a) and (2.6) we show 
the lowest eigenvalues versus <f> for {7 = and U = 5 t, respectively; k denotes the interplanar quasi-momentum 
quantum label. At r = the ground state energy is ^c2 3 o 12 (8) = ^Cu0 4 (^) F ^-Eq^q (2) an( ^ the ' ow ener gy 
sector derives mainly from the tensor product of the ground states of three independent CuCU's with 4, 2 and 
2 holes (which is the foundamental multiplet) and 3, 3 and 2 holes (which is the lowest lying excited multiplet 
separated by a gap Ac u o 4 (4)). For (7 = 0, see Fig. (2. a), there is no pairing in CUO4 and indeed the ground 
state energy is linear in the field at small fields (normal Zeeman effect). The lowest state is k = 2 throughout. 
Interestingly, CU3O12 concatenated with half a flux quantum would be diamagnetic, but the absence of a 
second minimum shows that it would be Larmor diamagnetism. By contrast, at [7=5 t, when pairing in CUO4 
is about optimum, see Fig.(l), the k = 2 state is lowest in the central sector, k = is the ground state at 
4> —> while k = 1 is lowest as <f) — > <f>o, see Fig.(2.&); this produces level crossings and the superconducting flux 
quantization; there is a central minimum when the system swallows a half quantum of flux while, as we verified, 
Acu 3 Oi 2 < 0. Remarkably, one also observes superconducting quantization of a magnetic flux orthogonal to 
the plane [12]. With increasing U/t, the binding energy of the pair starts decreasing and eventually vanishes 
for U ~ 34.77 t; we have found that at this point the flux quantization returns normal and the system behaves 
like a paramagnet. Even at optimal U, the side barriers are depressed by increasing r; at r ~ 0.1 1 only small 
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remnants remain; for still larger interplanar hoppings the superconductivity is removed and a pattern similar 
to Fig. (2. a) prevails. For larger |A|, a smaller supercurrent would be necessary to screen the half flux quantum 
and the critical r which kills the superconductivity should be expected to grow larger. 

In order to better analyze the results physically and extend them qualitatively to arbitrary graphs, we 
obtained an effective Hamiltonian by the cell-perturbation method with Ho, Eq.(3), the "cell-Hamiltonian" 
and H T , Eq.(4), the "inter-cell perturbation" and by taking into account only the low-energy singlet sector. 
We note that the cell-perturbation method was already used in Ref. [15] to support the original Anderson's 
conjecture [16] on the "low-energy equivalence" between the d — p model (proposed by Emery [17]) and the 
single-band Hubbard model. Despite the analogies with Ref. [15] (like the same cell-Hamiltonian and weak 
0-0 links between different cells) our inter-cell perturbation is different and, more important, it is the low- 
energy sector which differs (one needs to consider CUO4 units with 2, 3 and 4 holes to get bound pairs, in 
contrast with 0, 1 and 2 holes of Ref. [15]). 

For a general graph A, with 2|A| + 2p holes, we treat H T by a simplified second-order degenerate perturbation 
theory, since H T is a one-body operator. Each degenerate unperturbed ground state |4>o ) mav be labelled by 
a set S C A of units occupied by four holes; |<S| = p. The secular problem yields the eigenvalue equation 

&$\HT\*q)(* q \HT\&o )as> =e«s (8) 



A Cu o 4 (4) t-^t-f 

q S' 

where the sum has been truncated to the low-energy excited eigenstates involving Cu04 units with 2 < n < 4 
holes, all taken in their ground states |\I/q (a)), a = 1 . . . |A|. The amplitude as = a(«i, . . . , a p ) is totally 
symmetric with respect to permutations of the distinct indices ai, . . . , a p . Letting C(a) = {/3 £ A : r a /3 ^ 0}, 
after some algebra Eq.(8) may be written in the form: 



ea(ai, . . . , a p ) = >J 

3=1 0eC(a,) 



Bose Bosc -|— r 

Tp,a 3 o(ai, ..,aj-i,/3,a j+ i, ..,a p ) - \T 0:Ol . ]_]_(! - ^ a Ja(ai, 



i^3 



(9) 



This is a Schrodinger equation for p hard-core bosons hopping with an effective hopping integral T a a 
(r^) 2 /A Cu o 4 (4), with 



rf B = (*< 2) WI ® (¥ 4) (f3)\H T \¥ 3) (a)) ® |^ 3) (/3)>. (10) 
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In Fig. (3. a) we show the trend of (l^a^ |/|tq:/3 |) versus U/t; we note that the ratio decreases monotonically. 
In Eq.(9), the first term in the r.h.s. describes hole pair propagation, e.g. from unit ay to an unoccupied 
unit (5; in the second sum, the system gets back to the initial state after virtually exploring unit (5; the term 
— Sp ai ) takes into account that if f3 is one of the occupied units, the particle cannot move toward it. 

Due to the minus sign, the term in |7a iaj | represents pair-pair repulsion. 




FIG. 3.: (a) (|r^|/|r a(3 |) 2 versus U/t. (b) Results of Eq.(9) for Cu 3 Oi 2 with r = 0.001 t, U = 5 t. Lowest-energy 
eigenvalues labeled by their interplanar quasi-momentum are shown versus flux <fi. All energies are in t units. 

In Fig.(3.&) the superconducting flux-quantization for the |A| = 3 ring is reported as reproduced by solving 
Eq.(9); it agrees well both qualitatively and quantitatively with the numerical results of Fig.(2.b), thus con- 
firming the above approximation. More data and a fuller account of the low-energy theory will be presented 
elsewhere. 

In conclusion, we used a set of Cu04 units connected by weak 0-0 links to model interplanar coupling 
and c-axis superconductivity in Cuprates. The results show that the system with two holes in each unit is a 
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background such that inserting 2p holes one gets p pairs, bound by the repulsive interaction. The bound pairs 
propagation is well described by Eq.(9). We found analytically the superconducting flux quantization in the 
ring-shaped systems and confirmed this finding numerically for the 3-unit ring (f ,863,225 configurations). To 
this end, we introduced a novel exact-diagonalization technique, which reduces the size of the matrices that 
must be handled to the square root of the overall size of the Hilbert space. Actually, real systems contain 
also vertical links via the orbitals of the apical oxygens. We expect that the inclusion of these hoppings does 
not change qualitatively the results since they do not contribute to the propagation of the bound pair in the 
lowest-order approximation. 
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